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The nature of the chiral phase transition in lattice QCD is studied for the 
cases of 2, 3 and 6 flavors with degenerate Wilson quarks, mainly on a lattice 
with the temporal direction extension N t = 4. We find that the chiral phase 
transition is continuous for the case of 2 flavors, while it is of first order for 
3 and 6 flavors. 



1 Introduction 



In this article we investigate the nature of the chiral phase transition in lat- 
tice QCD for various numbers of flavors Np in the case of degenerate quarks, 
taking the Wilson formalism of fermions on the lattice. The determination of 
the order of the chiral transition, in particular for Np = 2, is a first step to- 
ward the understanding of what happens in the QCD phase transition, which 
is supposed to occur at the early Universe and possibly at heavy ion colli- 
sions. It is also important to compare the results for various number of flavors 
with theoretical predictions. Because the Wilson formalism of fermions on 
the lattice is the only known formalism which possesses a local action for 
any number of flavors, it is urgent to investigate the finite temperature chiral 
phase transition for various number of flavors with Wilson quarks. 

For each Np there are two parameters: ft (= 6/g 2 ; g is the gauge cou- 
pling constant) and the hopping parameter K. Because chiral symmetry is 
explicitly broken in the Wilson formalism even for a vanishing bare quark 
mass, we have to begin with the identification of the chiral limit. We define 
the quark mass through the axial Ward identity [I], 0] 

2m q < | P 1 7r >= — run < | A4 1 n > 

where P is the pseudoscalar density and A4 the fourth component of the local 
axial vector current. At zero temperature the pion mass vanishes at the point 
which is almost identical to that where the quark mass vanishesp], §]. It has 
been further shown that the value of the quark mass does not depend on 
whether the system is in the high or the low temperature phase through 
simulations at ft — 5.85 in the quenched QCD || and at ft — 5.5 for the 
Np = 2 case||. Therefore we can identify the chiral limit by the vanishing 
point of the quark mass m q . This defines a curve K c (ft) in the ft - K plane. 
Alternatively we may use the vanishing point of m 2 in the confining phase. 

At small ft region (ft < 5.0) where we have made simulations in this 
work, m q in the deconfining phase does not agree with that in the confining 
phase and the proportionality between the m q in the deconfining phase and 
m 2 in the confining phase is lost, contrary to the situation for ft > 5.5 
mentioned above. We interpret this as a lattice artifact. On the other hand, 
the proportionality between the m q in the confining and m 2 in the confining 
phase is well satisfied. Therefore we adopt as K c (ft) the vanishing point of m q 
in the confining phase K c (m q ) or that of m 2 in the confining phase K c (mfy. 
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The values of K c (f3) show a slight dependence on the choice of K c (m q ) or 
K c (m%)- They also slightly depend on Np and N t . Those at (3 = 4.5 are 
listed in Table [TJ. We find that the difference due to N F is of the same order 
of magnitude as that due to the difference of definition. The K c (m%) for 
N F = 2 for various /3's are listed in Table |[ 

The temperature on a lattice is given by 1/N t a, where N t is the lattice 
size in the temporal direction and a is the lattice spacing. The location of 
the transition or crossover K t (f3) from the high temperature regime to the 
low temperature regime, for a given (3 and a fixed temporal lattice size N t , is 
identified by a sudden change of physical observables such as the plaquette, 
the Polyakov line and screening hadron spectrum. 

The chiral transition occurs at the crossing point of the K t ((3) and K c (/3) 
line. However, whether the transition K t line with a fixed N crosses the 
K c line is not trivial. As was first noticed by Fukugita et a/.||, the K t line 
creeps deep into the strong coupling region. Therefore before discussing the 
order of the transition, we have to address ourselves to the problem whether 
the chiral limit of the finite temperature transition exists at all. 

In a previous paper [0] we showed that when Np > 7, the K t line does 
not cross the K c line at finite (3. In this article we identify the crossing point 
for the case Np < 6 and determine the order of the transition. Preliminary 
reports are given in M. We take the strategy of performing simulations on 
the line K c (f3) starting from a j3 in the high temperature phase and reducing 
(3. We call this method u on-K c " simulations. The number of iteration N inv 
needed for the quark matrix inversion, in general, provides a good indicator 
to discriminate the high temperature phase from the low temperature phase. 
The use of A^" inv as an indicator is extremely useful u on-K c " , because A^ inv is 
enormously large on the K c line in the confining phase, while it is of order 
several hundreds in the deconfining phase. This difference is due to the fact 
that there are zero modes around K c in the low temperature phase, while 
none exists in the high temperature phase 0, |^]. Therefore there is a sudden 
drastic change of A^ inv across the boundary of the two phases. Combining this 
with measurements of the Polyakov loop, the plaquette and hadron screening 
masses, we identify the crossing point (3 ct of the K t (f3) line with the K c (f3) 
line. We also check that the crossing point thus determined is consistent with 
an extrapolation of the line K t (f3) toward the chiral limit. From the behavior 
of physical quantities toward (3 ct) we are able to determine the order of the 
chiral transition. 
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2 Simulations 

We mainly perform simulations for the case of N t = 4 with the spatial size 
8 2 x 10. To study the Nt dependence for the Np = 2 case, we also make 
simulations at Nt = 6 and at N t = 18 with 12 3 and 18 2 x 24 spatial lat- 
tices, respectively. We generate gauge configurations by Hybrid Monte Cairo 
(HMC) algorithm for Np = 2 with a molecular dynamics step chosen in such 
a way that the acceptance rate is about 80 - 90%, while for Np > 3 by 
the hybrid R algorithm with the molecular dynamics step At = 0.01, unless 
otherwise stated. The inversion of the quark matrix is made by a minimal 
residual method or a conjugate gradient (CG) method. When the hadron 
spectrum is calculated, the lattice is duplicated in the direction of lattice 
size 10 or 12. We use an anti-periodic boundary condition for quarks in the 
t direction and periodic boundary conditions otherwise. The statistics is in 
general total r = several hundreds, and the plaquette and the Polyakov loop 
are measured every simulation time unit and hadron spectrum is calculated 
every 5t = 10 (or less depending on the total statistics). When the value of 
(3 is small, the fluctuation of physical quantities are small and therefore 
we think the statistics is sufficient for the purpose of this article. 

3 Results 

3.1 N F = 2 

Now let us discuss the results for the Np = 2 case at N t = 4. The _fQ(/3)'s 
obtained by various groups |l(J |[ are plotted in Fig. [I] together with the K c {(3) 
line[[Ll[]. In order to confirm the existence of the crossing point we take the 
largest (farthest) K c , that is K c (m^) for Np = 2 in Table |] and interpolated 
ones for u on-K c " simulations. We first perform u on-K c " simulations by the 
R algorithm to identify the crossing point, because it is very time consuming 
due to low acceptance rate for the HMC algorithm on K c in the confining 
phase. We find that when (3 > 4.0, A^ inv stays around several hundreds, 
while (3 < 3.9 it increases with r and exceeds several thousands (see Fig. |2|). 
Therefore we identify the crossing point at (3 ct ~ 3.9 - 4.0. This (3 ct is 
consistent with a linear extrapolation of the K t line as is shown in Fig. [TJ 
Our results for K t is summarized in Table ^. 
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Then we repeat u on-K c n simulations by the HMC algorithm for (3 > 4.0 
in order to measure physical observables. A^ inv for (3 > 4.0 plotted in Fig. |] 
is for the HMC algorithm, which is similar to that for the R algorithm. The 
At toward f3 ct should be taken small in order to keep the acceptance rate 
reasonably high (for (3 = 4.0, 4.1 and 4.2 we use At = 0.002, 0.005 and 0.005 
to get acceptance rates 0.91, 0.79 and 0.93, respectively). The value of m 2 
thus obtained decreases toward zero as the chiral transition is approached 
(see Fig. D. 

We find no two-state signals around (3^. This is in sharp contrast with the 
Nf = 3 and 6 cases where we find clear two-state signals at (3 ctl as discussed 
below. This together with the vanishing m 2 toward (3 c t indicate that the 
chiral phase transition is continuous for Nf = 2. 

The results by "on-i^ c " simulations on the N t = 6 lattice are similar to 
those on the N t = 4 lattice. The estimated transition point is (3 ct ~ 4.0 
- 4.2. The value of m 2 plotted in Fig. |3] again decreases toward zero as (3 
approaches /3 c t. For N t = 18 with the spatial size 18 2 x 24, we find that the 
transition is at (3 ct ~ 4.5 - 5.0. Although the spatial size is not large enough, 
this result suggests that the shift of (3 ct with Nt is very slow. 

3.2 N F = 3 

Now we perform u on-K c " simulations for Np = 3. In order to confirm the 
existence of the crossing point we take the largest (farthest) K c , that is 
K c (m^) for Nf = 2 at /3's we have studied, since this is the most stringent 
condition for the existence of (3 ct . We use them and interpolated ones here, 
and also for Np = 6 interpolated ones with K c = 0.25 at (3 — 0. 

Fig. |] shows iVi nv as a function of the molecular-dynamics time r for 
several values of /3's. When (3 > 3.1, A^" inv is of order of several hundreds, 
while when (3 < 2.9, iVj nv shows a rapid increase with r. At (3 = 3.0 we see 
a clear two-state signal depending on the initial condition: For a hot start, 
iVi nv is quite stable around ~ 800 and m 2 is large (~ 1.0). On the other 
hand, for a mix start, A^ nv shows a rapid increase with r and exceeds 2,000 
in t ~ 20, and in accord with this, m 2 decreases with r. 

The value of m 2 is plotted in Fig. At (3 = 3.0 we have two values for 
m 2 depending on the initial configuration. The larger one obtained for the 
hot start is of order 1.0, which is a smooth extrapolation of the values at 
/3 ~ 3.1 - 3.2. The smaller one is an upper bound for m 2 for the mix start. 
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We note that the result of (3 c t ~ 3.0 is consistent with an extrapolation 
of K t points listed in Table |3] as is shown in Fig. [l|. Thus we identify the 
crossing point at (3 ct ~ 3.0(1). With the clear two-state signal we conclude 
that the chiral transition is of first order for Np = 3. 

3.3 N F = 6 

Finally let us discuss the results for Np = 6. Our previous study [0] implies 
that this number of flavor is critical for the existence of the crossing point 
and therefore it is important to establish the existence of the crossing in 
this case. Overall features of the transition for Np = 6 are very similar to 
those for Np = 3 except for the position of /3 c t, which moves to a smaller 
(3 as expected. Fig. | shows that iVi nv "on-i^ c " stays at several hundreds 
for (3 > 0.4 and for a hot start at (3 = 0.3. On the other hand, A^ inv grows 
rapidly with r and exceeds 5,000 for (3 < 0.2 and for a mix start at (3 = 0.3. 
In accord with this, we have two values of m\ at (3 = 0.3 (Fig. 0). Therefore 
we identify the crossing point at (3 ct ~ 0.3(1) and conclude that the chiral 
transition is of first order for Np = 6. This (3 c t is consistent with a linear 
extrapolation of the K t line (K t = 0.245 - 0.2475, 0.235 - 0.237, and 0.166 - 
0.168 for (3 = 0.5, 1.0, and 4.5, respectively). 



4 Conclusions 

It should be noted that our main results obtained in this work that the 
chiral transitions for Np = 3 and 6 are of first order, while it is continuous 



for Np = 2, are consistent with the prediction based on universality [12]. Our 
results with Wilson fermions are also consistent with those with staggered 
fermions|l3]. 



The simulations on the N t = 4 and 6 lattices and the N t = 18 lattice 
have been performed with HITAC S820/80 at KEK and with QCDPAX at 
the University of Tsukuba, respectively. We would like to thank members 
of KEK for their hospitality and strong support and the other members of 
QCDPAX collaboration for their help. This project is in part supported by 
the Grant-in- Aid of Ministry of Education, Science and Culture (No. 62060001 
and No.02402003). 
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N t 


= 4 


N t = 8 


Np 




K c (m q ) 


K c {ml) K c (m q ) 


2 


0.214(1) 


0.210(1) 


0.212(1) 0.209(1) 


3 


0.209(1) 


0.204(1) 




6 


0.205(2) 


0.200(1) 





Table 1: Chiral limit K c at /3 = 4.5 determined by vanishing m 2 and m q . 
Values of m 2 and m q in the confinig phase are linearly extrapolated in 1/K 
using data from K — 0.16 - 0.18 for Np = 2 and 3 and K — 0.15 - 0.165 for 
Np = 6 (because K t = 0.167(1) for N F = 6 at N t = 4). The spatial lattice 
size is 8 2 x 10. 



p 


3.0 


3.5 


4.0 


4.3 


K c (ml) 


0.235(1) 


0.230(1) 


0.223(1) 


0.218(1) 


K c (m q ) 


0.230(1) 


0.226(1) 


0.218(4) 


0.214(1) 



Table 2: K c for A^^ = 2 determined on an 8 2 x 10 x 4 lattice. The results for 
(3 = 4.5 are given in Table |I|. 





Np = 2 






Np = 3 


p 


K t 


P 


K t 


4.3 


0.207- 


0.210 


3.0 


> 0.230 


4.5 


0.200 - 


0.202 


4.0 


0.200 - 0.205 


5.0 


0.170 - 


0.180 


4.5 


0.186 - 0.189 


5.25 


0.160 - 


0.165 


4.5* 


0.186 - 0.189 








4.7* 


0.179 - 0.180 








5.0 


0.166 - 0.1665 








5.5 


0.1275 - 0.130 



Table 3: Finite temperature transition/crossover K t for Np = 2 and 3 ob- 
tained on an 8 2 x 10 x 4 lattice (data with * obtained on a 12 3 x 4 lattice). 
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Figure 1: Phase diagram for iVp = 2 and 3. Filled symbols are for K c (m\) 
and K c (m q ). Open symbols are for K t (N t = 4) for Np = 2 while open circles 
with cross for Np = 3. Circles are our data. Lines are for guiding eyes. 
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Figure 2: Molecular-dynamics time history of A^ inv for N F = 2 on K c 
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Figure 3: for N F = 2 on K c . 
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Figure 5: for N F = 3 on K c . 
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Figure 6: Time history of N inv for N F = 6 on K c . 
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Figure 7: ml for N F = 6 on K c . 
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